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ABSTRACT 

Context. Circumstellar debris disks older than a few Myr should be largely devoid of primordial gas remaining from the protoplane- 
tary disk phase. Tracing the origin of observed atomic gas in Keplerian rotation in the edge-on debris disk surrounding the ~ 12 Myr 
old star y3 Pictoris requires more detailed information about its spatial distribution than has previously been acquired by limited slit 
spectroscopy. Especially indications of asymmetries and presence of Can gas at high disk latitudes call for additional investigation to 
exclude or confirm its connection to observed dust structures or suggested cometary bodies on inclined eccentric orbits. 
Aims. We set out to recover a complete image of the Fei and Can gas emission around /3Pic by spatially resolved, high-resolution 
spectroscopic observations to better understand the morphology and origin of the gaseous disk component. 

Methods. The multiple fiber facility FLAMES/GIRAFFE at the Very Large Telescope (VLT), with the large integral-field-unit 
ARGUS, was used to obtain spatially resolved optical spectra (from 385.9 to 404.8 nm) in four regions covering the northeast and 
southwest side of the disk. Emission lines from Fei (at 386.0 nm) and Can (at 393.4 and 396.8 nm) were mapped and could be used 
to fit a parametric function for the disk gas distribution, using a gas-ionisation code for gas-poor debris disks. 

Results. Both Fei and Can emission are clearly detected, with the former dominating along the disk midplane, and the latter reveal- 
ing vertically more extended gas. The surface intensity of the Fe i emission is lower but more extended in the northeast (reaching 
the 210 AU limit of our observations) than in the southwest, while Can shows the opposite asymmetry. The modelled Fe gas disk 
profile shows a linear increase in scale height with radius, and a vertical profile that suggests dynamical interaction with the dust. We 
also qualitatively demonstrate that the Can emission profile can be explained by optical thickness in the disk midplane, and does not 
require Ca to be spatially separated from Fe. 
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1. Introduction 

The transformation of dusty and gaseous disks around young 
stars into planetary systems involves a gradual clearing of 
the disk by several processes. Dense and gas-rich protoplane- 
tary envelopes are affected by strong stellar radiation ( Gorti & 
|HoUenbach||2009] l and winds ([Lovelace et al.|2008), that co m- 
pete with viscous accretion (Lynden-Bell & Pringle] [l974| l to 
open up growing inner disk holes, and cause complete pho- 
toevaporation of the gaseous disk component within -10 Myr 
(jZuckerman et al.[T9 95 ; Hais ch et al.|200l| |Jayaw ardhana et al. 
2006;; |Pascucci et al.|2006[ ). This also sets the time frame for gi- 
ant planet formation, because it restricts the amount of gas avail- 
able to be swept up by growing planetary embryos. Primordial 
dust grains are blown out from the system or spiral onto the 
star due to Poynting-Robertson (PR) drag and are cleared from 
the disk within a few thousand years. The left-over debris disk 
is defined by a low (<10%) gas-to-dust ratio, is mostly opti- 
cally thin, and often observable as excess mid- and far-infrared 
(IR) emission in the spectral energy distribution (SED), origi- 
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nating from heated dust grains that are being continuously pro- 
duced in collisions between planetesimals (Dominik & Decin 



2003 Wyatt et al. 2007 1. Decrease in this steady-state colli- 



sional processing and continous removal of submicron sized 
dust by radiation pressure leads to a decline in observed IR 
excess with time. Most of the dust in a standard debris disk 
is believed to be composed of porous silicate grains contain- 
ing Si, Mg, and Fe, likely combined with oxygen into forsterite 
(Mg 2 Si0 4 ), enstatite (MgSi0 3 ), olivine ([Mg,Fe] 2 Si0 4 ), and py- 
roxenes ([Mg,Fe]SiC>3), as found in chondritic meteorites (see 
review by Henning 2010). Beyond the ice condensation (snow) 
line, grains would contain a cometary fraction of H2O, CO2, CO, 
NH3, CH4, and N2 ices, although there is a possibility that photo- 
sputtering might strip grains from their ices beyond the snowline 
(Grigoriev a et al.|2007 1. 



Observation of gas in such evolved non-primordial disks 
is not expected a priori. It is of course routinely observed in 
younger systems in the protoplanetary-disk stage. As an exam- 
ple, abundant molecular species (e.g. CO and H2) in young cir- 



cumstellar disks have been detected in both emissi on ([Najita 
l et al.|2003||Herczeg et al.|2006) and absorption ( |Redfield|2007[ 
Roberge & Weinberger 2008 ). However, in most disk environ- 
ments, especially around early-type stars, molecules like H2O, 
OH, and CO are not expected to survive, instead producing an 
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abundance of atomic oxygen and carbon by photodissociation 
(e.g. Kamp & Bertoldi 2000 1. Elements with strong resonant 
transitions in the ultraviolet (UV) and optical, especially metals, 
are then easily blown out from the system by radiation pressure. 
Thus, the detection of substantial quantities of gas in older sys- 
tems that have presumably made the transition to debris disks, 
like the 30 Myr-old HP 21997 (poor et al.|201 1) and 200 Myr- 
old <x Her ( |Chen & Jura 2003]!, is surprising and points to un- 
known gas retention or production mechanisms. If gas is being 
produced, it can give us information on the composition of solid 
bodies in the disk, and even small amounts of gas can affect the 
dynamics of dust grains (Thebault & Augereau 2003J |Krivov| 



et al. ||2009| l . j3 Pictoris is the only currently known debris disk 
that contains gas and is close enough to be spatially resolved 
at optical and IR wavelengths, allowing detection of resonantly 
scattered light and thus a more detailed study of gas composition 
and distribution. 

The first indication of circumstellar material around the 
~12Myr A6V star /?Pic came from detection of thermal IR 
excess by the Infrared Astronomical Satellite (IRAS, |Aumann| 
1985) 1, and subsequent ground-based coronagraphic imaging re- 
vealed a flared, nearly edge-on dust disk in scattered optical and 
near-IR light ( [Smith & Terrile|1984|[Paresce & Burrows|1987| l, 
extending out to 1800AU. Following its classification as a de- 
bris disk star (or so-called Vega-like star, from the prototype 
IR excess star Vega, |Backman & Paresce|1 993 ), (3 Pic has been 
extensively studied and the minimum mass of cold dust con- 
tained within its disk determined from submillimeter (submm) 
and millimeter observations to 3-10 lunar masses (|Zuckerman 
& Becklin|19"93HHolland et al || 19981 [Dent et al.||2000[ |Liseau| 
et al.|2003[ [Nils son et al .|2009|l. In ad dition to the collisionally 



produced dust (|Backman & Paresce 1993 Artymowicz 1997 



Lagrange et al.||2000[ Zuckerman 2001 1, several lines of evi 
dence suggest that evaporating comets on eccentric orbits, so 
called falling evaporating bodies, FEBs, are supplying dust to 
the jSPic debris disk ( |Lecavelier Pes Etangs et al.|1996l|1998 1. 
One comes from spectroscopic signatures of intermittent red- 
shifted absorption lines ( Lagrange -Henri et al.|19 88 1 Beust et a 



1990 |Vid al-Madjar et al. 1994 and references therein) and a 
second from emission by crystalline silicates (at 9.7, 28, and 
33.5 //m) and olivines (at 11. 3 /mi) (|Telesco & Knacke|[l991 



Knacke et al.p993| |Aitken et al.|1993||Wei nberger et al. 2003 
Okamoto et al.||2004[ |Chen et al.||2007) , similar to that ob 
served in e.g. comet Halley. This points to a cometary-like 
grain composition of either crystalline olivine (~55%), pyroxene 
(~35%), and other silicates, or a mix of 95% amorphous and 5% 
crystalline olivine ( |Artymowicz 1997| and references therein). 
Dissimilarities between disk features in IR emission and scat- 
tered light images of the disk imply that two distinctly differ- 
ent grain populations contribute, perhaps low-albedo refractory 
grains (like the hot, > 300 K, silicates that are unmistakably re- 
sponsible for the 10 pm emission) and high-albedo icy grains. 

Information about the spatial distribution of circumstellar 
gas is more limited than the dust, since it is more difficult to 
observe. Several species of atomic gas have been detected in ab- 
sorption, due to the disk's favorable edge-on orientation. The 
strong circumstellar CanK absorption line, e.g., was detected 
early, together with weaker NaiD lines ( |Vidal-Madjar et aL 
1986[ Hobbs et al.|1 985 ), and several neutral and singly ionised 
metallic species have been observed in both optical and UV ab- 
sorption spectroscopy (Lagra nge et al.|1 998). Surprisingly, these 
appear to be very stable (apart from the intermittent redshifted 
components) and do not show any radial velocity, which would 
be expected from acceleration by radiation pressure. 



Olofsson et al. (2001) performed high-resolution slit spec- 
troscopy along the disk of j3 Pic and made the first spatially re- 
solved detection of circumstellar atomic gas around any star. The 
NaiD emission line (doublet) at 5895 A was found to be reso- 
nantly scattered from 30 to 140 AU, with an average Na i column 
density of 10 15 cuT 2 and a velocity pattern clearly indicating gas 
in Keplerian rotation. To explain the absence of any detected 
radial motion - considering the fact that the force from radia- 
tion pressure felt by Na i throughout the ft Pic disk is 250 times 
stronger than the gravitational force - Olofsson et al. ( 2001| l pro- 
posed fast ionisation by the intense stellar UV radiation, which 
does not allow significant time for acceleration of neutral Nai 
before conversion to Na n ions (which lack strong transitions in 
the wavelength region dominating the stellar SED). This model 
was recently improved by Brande ker] ( |20 1 1 ) . In VLT/UVES ob- 
servations over a wider spectral range and with a spatial and 
spectral resolution twice as high, Brandeker et al. ( 2004[ ) de- 
tected 88 emission lines originating from Fei, Nai, Can, Nil, 
Nin, Tii, Tin, Cri, and Cm. Spatial information from the slit 
placed at four positions across the disk (at 60 and 120 AU, on 
either side) and four positions along the disk revealed a surpris- 
ingly extended and asymmetric gas distribution. The strong Fe i 
and Na i emission reaches along the disk from 13 out to 323 AU 
from the star, displaying a NE-SW asymmetry similar to that 
observed in the dust, with a gradual NE radial decline and a 
sharp cutoff in the brighter (out to 75 AU) SW emission at 150- 
200 AU. Also in accordance with dust observations, an inner 5° 



tilt was found by [Brandeker et aL (20041. The radial thickness 
increase of the gas disk appears to be significantly higher than 
that of the dust (almost twice the scale-height of the dust disk 
at 116AU). Can H- and K emission can be traced to heights 
of 77 AU above the mid-plane, but falls off closer to the mid- 
plane. As in the observations by Olofsson et al. (2001), the out- 
ward radial motion of gas is much slower than would be ex- 
pected from acceleration due to radiation pressure. This points 
to a braking agent that cannot be Hi and H2, which have ob- 
servational upper limits in column density inconsistent with ob- 
served velocities, or perhaps an unknown gas production mech- 
anism (Brande ker et al.||2004| l. [Fernandez et al.| < (2006] > studied 
gas braking and production processes, concluding that decelera- 
tion of ions through Coulomb interaction with other ions (self- 
braking) could be sufficient to brake the gas, but only if C in 
the disk gas is at least five times overabundant compared to so- 
lar elemental abundances. Spectroscopic observations with the 
Far Ultraviolet Spectroscopic Explorer (FUSE) (Roberge et al. 
2006 ) did find C 11 and C m measured in absorption to be signifi- 
cantly more abundant compared to elemental abundances of so- 
lar system material (viz. the Sun, carbonaceous chondrites, and 
dust from comet Halley), and could thus partly solve the radial 
velocity conundrum. More detailed modelling of the braking by 
Brandeker ( 201 l| l and observations by the PACS instrument at 



the Herschel Space Observatory (Brandeker et al. 2012) indi 
cated that the C abundance over metallic elements such as Na 
and Fe might be even higher, of about 400x solar abundance. 

Because of gas-removal mechanisms, the gas is probably not 
primordial (Fernan dez et al.| [2006), but replenished. Suggested 
mechanisms include (1) gas released from star-grazing comets 
( |Beust &~V aliron 2 007|l, (2) photo-desorption from circumstel- 
lar dust grains ( Chen et aT[ 2007 ), and (3) vaporisation of dust 
by collisions with high-velocity yS-meteoroids (Czechowski & 
Mann 2007). Since these different mechanisms should produce 
different spatial distributions of gas, one way to potentially dis- 
tinguish between them is to map the spatial distribution of gas, 
which is the aim of this study. One would, e.g., expect most of 
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the gas released by comets to come from the inner regions of 
the system; the photo-desorption release of gas would be ex- 
pected to be oc A(r)/r 2 , where r is the distance to the star and 
A(r) is the dust area density (since the UV radiation responsible 
for photo-desorption drops as r 2 ); while the dust-dust collision 
gas production profile would need more detailed modelling (as 
shown by |Czechowski & Mann|2007| l. Once produced, the gas 



Observed positions and image scale of ARGUS array along |3 Pic disk 



may get distributed/removed, which is another effect that needs 
to be taken into account to properly predict the current spatial 
distribution of gas. This detailed modelling is beyond the scope 
of the present paper, however, which focuses on the observed gas 
profile. 

To obtain a complete map of the Fe i (at 386.0 nm) and Can 
(at 393.4 and 396.8 nm) emission throughout the disk of /?Pic, 
we used high-resolution integral field spectroscopy to cover the 
disk (Sects. [2] and [3J. The observed emission from Fei is mod- 
elled to obtain the spatial number density distribution of gas, 
combining a de-projection technique (Sect. |4.1[ ) with the ion- 
isation and thermal equilibrium code Ontario (Sect |4.2[ ). The 
obtained gas distribution is then used to compute the expected 
shape of the Ca n H emission profile, using a simple radiative 
transfer code to handle the optically thick H line (Sect, 
conclusions are summarised in Sect. [6] 



12 



10- 



8- 



4- 



0- 



-2- 



-6 



-10- 



-12 



5.1 1. Our 



0.52"/pixel resolution 
0.30"/pixel resolution 




4 2 -2 -4 
R.A. offset (arcsec) 



2. Observations and data reduction 

2.1. Observations with FLAMES/GIRAFFE at VLT 

The multi-object fiber facility FLAMES at the 8.2-m diameter 
Kueyen (UT2) telescope of the Very Large Telescope (VLT) ar- 
ray was used to obtain spatially resolved, high spectral reso- 
lution optical data of the ySPic disk, employing its GIRAFFE 
spectrograph with the ARGUS fiber configuration. ARGUS is 
a large integral-field-unit (IFU) composed of a rectangular ar- 
ray of 22x14 microlenses, sampling either 0.52"/microlens 
(with a total sky coverage of 11.5"x7.3") or 0.30"/nricrolens 
(6.6"x4.2"), and feeding the spectrograph through individ- 
ual fibers. The high-resolution (HR) dispersion grating (with 
/?=31300 or dAx 0.0126 nm) offers a total spectral range of 
370-950 nm accessed through various setups. We used the bluest 
HR setting with the central wavelength of the GIRAFFE grat- 
ing at 395.8 nm. Observations were carried out in service mode, 
spread out over a period ranging from November 2008 to 
January 2010 (due to stringent < 0.5" seeing requirements), with 
an average airmass of 1.17. The ARGUS array was placed at four 
positions orthogonally along the disk (two on each side of the 
star), partially overlapping the smaller scale image of the inner 
disk with the larger scale image of the outer disk (as displayed 
in Fig. [T| to check for consistency. At each position, exposures 
were taken at two position angles (rotated 180°) to reduce the 
risk of having any sky areas repeatedly falling on dead pixels. 
Table [T]presents observing dates, field centre positions, position 
angles (PA), pixel scales, and exposure times of the observations. 

2.2. Reducing FLAMES/GIRAFFE data 

The raw images containing dispersed spectra from all fibers (in- 
cluding sky- and simultaneous calibration exposures) were re- 
duced using ESO's recipe execution tool EsoRex with supplied 
GIRAFFE pipeline recipes. Calibration products (master bias, 
master flat-field with fiber localisation data, dispersion solution 
and line identification from simultaneous ThAr arc lamp spec- 
tra, and instrument response from standard star observations) 
were created and then applied in the reduction of science frames 



Fig. 1: Scale and position in the sky of the ARGUS array during 
the observations. The northeast (NE) and southwest (SW) side 
of the optical disk plane (marked with a black line) was cov- 
ered, omitting the area closest to the star, where the stellar PSF 
and scattered light would be strong, giving low signal-to-noise 
ratio. Cyan- and magenta-coloured areas represent the higher 
(0.30"/pixel) and lower (0.52"/pixel) spatial resolution settings 
used for the inner and outer disk, respectively. 

Table 1: Observing log from VLT/FLAMES/GIRAFFE obser- 
vations of jSPic. The optical coordinates of the star are R.A. = 
05h47m 17.088s and Dec. = -51°03'59.44" (J2000). 
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Notes. <o) ARGUS field centre position. 



to retrieve flux- and wavelength-corrected spectra in each pixel 
of the reconstructed image. The final spectra were extracted by 
summing up over all pixels in the fiber width and were subse- 
quently rebinned to a linear scale of 0.005 nm spectral elements 
(from the raw ~0.0046 nm/pixel non-linear scale), ranging from 
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Table 2: Observed emission lines. 



w Pl/V \A H-£X3889 




390 392 394 396 398 400 402 404 
Barycentric wavelength, A (nm) 



Fig. 2: Total observed near-UV spectra of inner disk, comprised 
mainly of a stellar PSF component, with additions from gas 
emission and light scattered off circumstellar dust grains. 



385.9 nm to 404.8 nm, with a measured average resolution (from 
ThAr lines) full- width at half-maximum (FWHM) of 0.01 1 1 nm. 
Line-flux errors were estimated from the noise in the region 
surrounding each line, but it should be noted that the absolute 
photometric accuracy of the instrument response at these blue 
wavelengths are not better than ~20%, judging from the stan- 
dard deviation of obtained instrument response curves (see |A.1| 
in Appendix [A]). 

The data cubes were further processed in matlab, where 
barycentric wavelength correction and airmass correction was 
applied, and images were combined. Most cosmic ray hits could 
be identified and eliminated by sigma clipping and vertical 
(spectral dimension) interpolation directly in the unprocessed 
images. Remaining spikes in the spectra of the inner disk could 
be removed by modified median filtering, replacing values in the 
contaminated spectrum with the mean of the two additional ex- 
posures obtained (after normalisation), but the long single expo- 
sures of the outer disk presented a bigger problem. In that region 
we had to replace identified spikes with values interpolated from 
neighbouring bins, taking care to correct for occasional adjacent 
or extraordinarily wide spikes by comparing with the rotated ex- 
posure. The fact that the unresolved emission lines we were map- 
ping had a width comparable to the average width of the cosmic 
ray spikes added to the complication. 

2.3. Isolating the emission lines 

To extract the flux from the Fe i and Ca n emission lines (Table|2]i 
we had to subtract the dominant stellar point-spread function 
(PSF) and dust-scattered light contribution to the total spectrum 
(Fig.01. 

A spectrum (with FWHM resolution ~0.004nm and sam- 
pling of 0.001 nm/pixel) of /3Pic showing no signs of FEBs, 
which had previously been obtained from the HARPS spec- 
trograph at the ESO 3.6-m telescope, was convolved with a 
Gaussian of FWHM=0.0111 nm (measured average resolution 
from arc spectra) and rebinned to 0.005 nm, and used as a ref- 
erence stellar spectrum, /*. This was then fitted to the observed 
spectrum, / tot , surrounding each emission line region and in each 
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Notes. Data retrieved from NIST (http://physics.nist.gov). 



pixel individually by unconstrained non-linear minimisation of 
the function 



ZdUkW) = /totW) - (M 2 + M + h)U(X) - k 4 



(1) 



essentially multiplying it with a second-degree polynomial 
(which was interpolated across the line), but with compensation 
for offset level. This baseline was then subtracted to recover the 
gas emission lines. For the inner disk observations of Can we 
substituted the HARPS spectrum with a stellar spectrum derived 
from the mean of eight off-disk pixels with strong stellar PSF 
contribution, since redshifted Ca n components from FEBs were 
visible in these epochs. The total line flux was found by integrat- 
ing the spectral energy distribution over the approximated line 
width at the central wavelength position of the Doppler shifted 
line (assuming Keplerian rotation of the gas disk). These steps 
involved some tweaking of parameters (e.g. regions for baseline- 
fit and interpolation, and fitting function) to minimise the effects 
from noise and certain spectral features. Specifically, the Can 
H- and K emission lines are superimposed onto the broad stellar 
and narrow circumstellar Ca n H- and K absorption lines, mak- 
ing the fit of the stellar spectrum very sensitive to the chosen 
fitting region, and consequently affecting the subtracted base- 
line. Although this overlap made the reduction more difficult, it 
also helped by suppressing the background PSF noise of the Ca n 
emission. Fig. 3] shows emission lines of Fe i and Can H and K 
in two selected pixels, together with fitted baselines and regions 
used in the extraction. The baseline fitting regions extend 100 
and 500 wavelength bins on either side of the emission lines for 
Fei and Can, respectively. The position of the line peaks co- 
incides with expected Doppler-shifts due to Keplerian rotation 
relative to the systemic velocity of ySPic. Apparent redshifted 
offsets from systemic velocity in the Ca n H- and K absorption 
is a clear indication of the presence of FEBs. 

The signal-to-noise ratio (S/N) of the emission lines was es- 
timated by dividing the baseline-subtracted peak flux with the 
standard deviation, <x, of the surrounding spectral region. 

3. Results 

3.1. Disk seen in gas emission 

Resulting maps of FeiA3860 and Can H- and K emission in 
the disk of fiPic are presented in Figs. |4}j6] where the (b) pan- 
els display the intensity scaled with projected distance to the 
star squared to bring out detail in the outer regions. The opti- 
cal dust-disk midplane is marked by a blue line, and the tilted 
4.5° secondary disk with a cyan-coloured line. Fei A3860 emis- 
sion is clearly visible along the disk midplane, while Can H- 
and K emission is strongest in higher latitude areas on the NE 
and S W side. Both exhibit a highly asymmetric distribution; Fe i 
being fainter but more extended in the NE than in the SW, and 
Ca ii showing the opposite behaviour, in addition to asymmetries 
relative to the disk midplane. Fe i is detected at disk radii out 
to the 210 AU limit of our observations in the NE, while Can 
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39.1.65 



Barycentric wavelength. A (nm) Barycentric wavelength, A (ntn) 



(a) (b) 



Fig. 3: Example spectra of Fe i A3860 (a) and Ca n K (b), both sampled from the NE side of the disk. The regions used for baseline 
fitting and line integration are indicated, as is the systemic rest velocity of/3Pic (dashed vertical line). 




Projected distance from star (AU) Projected distance from star (AU) 



(a) (b) 

Fig. 4: Map of Fei emission, multiplied with the projected distance to the star squared in panel (b). The blue line indicates the disk 
midplane (defined to be at 30.75°) while the green line shows the orientation of the inner disk, as found in dust disk observations. 



appears to be reaching beyond the observed ~100AU latitudes 
on the SW side. Data from the innermost 2-3 columns (~20- 
30 AU) are very noisy due to strong contamination by light from 
the stellar PSF, therefore it is hard to tell if there is a decline 
in flux inward of the peaks at ~40 AU. However, inspection of 
the maps normalised with distance squared suggests that there 
is a peak in the column density of Fei at ~90 AU in the NE and 
~70 AU in the SW. Because Ca n H- and K emission is optically 
thick, we cannot immediately draw similar conclusions for Ca n, 
as discussed in detail in Sect. 15. II 

A brightness difference can be observed in the overlapping 
regions between observations of the outer and inner disks. This 
can be caused by a number of factors and is discussed in more 
detail in Appendix |A| 

Fig. [7] clearly shows how the emission contours from neutral 
Fe outlines the disk midplane, while singly ionised Ca emission 



appears weaker in that same region, instead dominating at high 
latitudes and increasing with radial distance, forming a charac- 
teristic butterfly shape. Note that this figure is only presented to 
highlight the relative positions of regions with dominating Fe i 
and Can gas emission, and that the contour levels of the inner 
and outer parts differ. 

4. Analysis 

Our primary interest is to find the spatial gas density distribution 
in the disk. The location of the gas gives clues to its origin, and 
can also put constraints on potential planet-disk interaction (as 
has indeed been done for the observed spatial dust distribution). 
Of particular interest is the Can observed at surprisingly high 
scale-heights above the disk midplane, and the possible mecha- 
nisms that may have placed the Ca there. 
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Fig. 6: As in Fig. g but for Ca n H. 



The problem of inverting observed flux into gas density is 
not trivial, because we are observing an essentially edge-on pro- 
jection of the disk emission, where the emission at each pro- 
jected distance is the contribution from many locations in the 
disk along the line of sight. If we furthermore consider the pos- 
sibility of an asymmetric disk, the problem of finding the gas 
density distribution from the observed projected flux becomes 
degenerate. To break the degeneracy, we assumed that the disk 
is cylindrically symmetric. Since we observe significant differ- 
ences between the NE and SW sides of the disk, we know that 
this assumption is of limited validity; we therefore made inde- 
pendent inversions of the two sides to get an impression of the 
uncertainties involved. To simplify the procedure even more, we 
divided the problem into two steps: 

1. De-project the Fei A3860 emission from the disk to obtain 
the spatial luminosity density as a function of cylindrical ra- 
dius r and height z (assuming cylindrical symmetry). 



2. Use the derived 2D luminosity density profile to derive the 
Fe gas density, with the help of the ionisation, thermal bal- 
ance, and level population statistical equilibrium code On- 
tario. 

The inversion from the Fe i A3860 luminosity density to the 
number density of Fe is complicated by the complex energy level 
structure of Fei. The level population is dominated by photo- 
excitation, hence the need for a non-LTE (local thermal equi- 
librium) solver (part of Ontario). Furthermore, the Fei A3860 
emission i s expected to be opti cally thin; for a line width of 
1 .5 km s _1 ( Crawford et al. 1994 1, the mid-plane column density 



of Fe i towards the star is 7X10 11 cm 2 ( Zagorovsky et al. 



2010) 



and conservatively assuming all atoms to be in the ground state, 
the optical depth of the A3860 transition becomes t ~ 0.1. 

In contrast, assuming a solar abundance ratio of Fe/Ca 
(Zagorovs ky et al.|2010 l, the Can H- and K lines are strongly 
optically thick (t ~ 100 in the mid plane). We therefore did 
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Fig. 7: Contour map of Fe i and Can H- and K emission. The contour levels of inner and outer panels differ and are just meant to 
give a clearer picture of the relative positions of dominant Fe i and Ca n emission regions. 



not attempt to use the Ca n lines to derive the gas density of Ca, 
since even small deviations from cylindrical symmetry in density 
will give large variations in emitted flux, making the inversion 
strongly degenerate. Instead, we used the spatial gas disttibution 
derived from Fe i, solved for the radiative ttansfer and computed 
the expected Ca n H- and K emission from the disk, in an attempt 
to check whether there is chemical segregation, as has been pre- 
viously suggested ( Beust & Valiron 2007 1 



Table 3: Best-fit parameters and errors. 



In the following subsections, the inversion steps are outlined 
in detail. 



4.1. Finding the Fe i \3860 luminosity density profile 

To find the de-projected luminosity density profile of Fe i A3860, 
we assumed cylindrical symmetry and fitted a few-parametric 
functional form to the observed flux distribution. In line with 
Brandeker et al. (j2004), we used a broken power-law for the 
radial dependence but inttoduced two new parameters for the 
height dependence (y) and the scale-height dependence (d>) on 
the midplane radius r. 





2 


1/2 






Lir,z) = k 


(r/r ) 2 » + (r/roF 


exp 


(h(r)) 


, (2) 



where a and /3 are the parameters describing the shape of the 
power-law on either side of the inflection point ro, and the verti- 
cal scale-height at each r defined as 



h(r) = ho 



''<> 



(3) 



The parameter 6 describes the 'flaring' of the disk. For <5 = 
we have a constant thickness disk, while for 6=1 the disk scale 
height increases linearly with radius. For 6 > 1 we have a flaring 
disk, i.e. the disk opening angle increases with radius. The pa- 
rameter y determines how quickly the density drops with height. 



Fe i A3860 luminosity 



Density 





NE 




SW 


NE 


SW 


Parameter 


value 


(T 


value 


a 


value 


value 


l (lO-'^ergs-'cm^) 


5.54 


1.10 


8.59 


1.60 






n (10 3 cm- 3 ) 










1.19 


1.76 


h (AU) 


16.4 


1.0 


16.9 


1.0 


21.8 


20.4 


r (AU) 


87.7 


4.9 


85.4 


3.3 


94.3 


83.9 


a 


4.98 


0.15 


6.88 


0.24 


3.00 


3.86 


P 


1.15 


0.28 


1.38 


0.29 


0.83 


0.93 


y 


1.50 


0.11 


1.57 


0.11 


1.26 


1.33 


5 


0.93 


0.07 


1.19 


0.10 


0.96 


1.21 



For a gas disk in thermostatic equilibrium we would expect the 
number density to drop as a Gaussian function, i.e. y — 2. The 
dust distribution, on the other hand, is best fitted by closer to an 
exponential (y w 0.5; Ahmic et al. 2009). 

For a given set of parameters, L(r,z) was projected onto a 
plane from edge-on to a resolution ten times better than the ob- 
servations. The projected flux distribution was then convolved 
with the seeing recorded for the observation by using a Gaussian 
of the same FWHM, and re-sampled to the grid of the obser- 
vations for comparison with observed data. By computing the 
residuals weighted by the noise, the best-fit values for the pa- 
rameters Iq, ho, ro, a,/3, y, and 6 were found by ^-minimisation. 
Fig.[8]displays maps of observed Fe i emission, estimated noise, 
best-fit degraded model, and weighted squared residuals of the 
model fit. Resulting parameter values from fitting of the gas-disk 
profile and associated errors are presented in Table [3] The errors 
were estimated using the Monte-Carlo method to make repeated 
fits to the data with added noise drawn from the calculated noise 
distribution in each pixel. 
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Fig. 8: Display of input data and results from ^-minimisation disk-modelling of the Fei emission. The rows of panels show, from 
the top down: observed Fe i emission maps, derived noise maps, degraded model maps, and maps of the weighted squared residuals. 
The reduced^ 2 values are 1.9 (NE) and 2.3 (SW). 



4.2. Finding the gas density from the luminosity density 

To find out how much gas we need to produce the inferred lu- 
minosity density of a Fe i A3860, we have to answer these ques- 
tions: 

1 . What is the number fraction of Fe in the gas? 

2. What is the fraction of Fe in the neutral, atomic state? 

3. What is the fraction of the Fe i that is radiating at A3860? 

Because f3 Pic is an A-type star, we expect its radiation field to 
dissociate almost all molecules in its disk, leaving only atomic 
and ionised gas components. This simplifies the estimate of free 
Fe atoms because we can ignore chemical interactions, and only 
need to consider the ionisation state. For this we used the Ontario 
code (Za gorovsky et al.]|2010 l, which self-consistently models 
the ionisation states and gas temperatures in debris disks around 



A- and F-type stars given a known stellar flux and gas/dust pro- 
file. It performs a full statistical equilibrium computation for 
several atomic species to find dominant absorption and emis- 
sion lines. The local gas temperature is determined by the bal- 
ance of heating and cooling mechanisms, and can, for the as- 
sumed optically thin gas, be treated in independent computation 
bins. The dominating heating processes for /? Pic are photoelec- 
tric heating by dust, and photoionisation of gas by stellar UV 
light, while cooling is dominated by fine-structure transitions (in 
particular from C n and O i, as recently detected in emission by 
the Herschel Space Observatory, Brandeker et al. 2012). A de- 
tailed description of the code can be found in Zagorovsky et al. 
( |20T0l >. 

By iterating Ontario with different gas densities, the density 
that best reproduced the given emission for each grid cell was 
found. The abundance of elements in the gas was assumed to 
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be solar ( |Grevesse & Sauval|[l998[ >, except for H (0.1% of so- 
lar), He (none), and C and O (400 times solar) as inferred from 
observations ( Brandeker et al.|2012 1. The resulting gas profiles 
(separately for the NE and SW sides) are presented in Fig. ^to- 
gether with a best analytical fit, using the same parameters as 
L(r, z) of Eq. [2] but with the luminosity density Zq replaced by a 
number density no, which is a parameter that refers to the number 
density of H had it been at solar abundance. The number density 
of any given element can then be found by multiplying «o with 
the abundance relative to H (had H been at solar abundance). For 
instance, the number density of Fe is given by npe = 2.8 x 10~ 5 «o 
and C by n c = 9.8 x 10~ 2 n . 



5. Discussion 

Looking at the spatial distribution, we confirm the strong asym- 
metry between the NE and SW sides of the disk, where the SW 
gas density seemingly drops much quicker with radius at a rate 
oc r~ 3 9 compared to the NE side oc r~ 30 . Neither side shows 
evidence for flaring, as 6 ~ 1. This means that the h/r ratio is 
fairly constant and close to ho/ro ~ 0.2. The height dependence 
with y ~ 1.3 is closer to exponential (y = 1) than the Gaussian 
(y = 2.0) one would expect from hydrostatic equilibrium. This 
could be an indication that the gas dynamics is significantly af- 
fected by the dust grains, as found by Fernandez et al. (2006]! to 
be expected for an ionised gas interacting with charged dust. 



5. 1 . The Ca n profile 

In contrast to the detected Fe i emission (Fig. [4), the Ca n H- and 
K emission (Figs.[5Jj6]) is observed far away from the midplane. 
An interesting question is whether there is a mechanism prefer- 
entially transporting Ca to far heights above the midplane (as ar- 
gued by e.g. Beust & Valiron|2007| l, or if there is another reason 
for Ca ii emission to be preferentially emitted there. A significant 
difference between Ca n and Fe i (and the other spatially mapped 
species Na i from Brand eker et al.|2004[ ) is that, since the metals 
are all strongly ionised, the number density of Can is 15-500 
times higher than the neutral Fe i (and Na i). Furthermore, the H- 
and K transition probabilities of Can are lOx higher. Together, 
this means that while the neutral transitions are optically thin, the 
Ca ii lines are optically thick. A reason for the paucity of mid- 
plane emission could therefore be that the star is obscured by the 
gas, and that the only photons being scattered from the H- and K 
transitions in the midplane either come from secondary sources 
(e.g. photons scattered down from higher altitudes) or from far 
out in the wings, where the lines are still optically thin. 

To test the idea that the optical properties are responsible for 
the difference in appearance between the Fe i A3860 and Ca n H- 
and K emission maps, we implemented a simple radiative trans- 
fer code to compute the emerging emission map. Because we 
expect the lines to be optically thick, the line shape is important 
as the optical depth will vary over the line profile. Furthermore, 
because the gas is orbiting the star, there will be a velocity field 
that will Doppler-shift the line profile depending on the velocity 
projected on the line of sight. In detail, our model is based on 
the following assumptions and parameters: 

1. The spatial distribution of gas is assumed to be the one in 
Table 3] inferred from the Fei A3860 emission (we assume 
different profiles for the NE and SW sides). The number den- 
sity of Can is then computed using Ontario with the same 
parameters as used when inverting the Fe i emission profile . 



2. For the velocity field, circular Keplerian orbits are assumed 
with a dynamical mass equal to Md yn = 1.4 M® (Olofsson 
|et al.|2001) . 

3. As an intrinsic line profile, a Voigt profile is assumed with 
broadening parameter b = 1.5kms _1 . This is close to, but 
slightly below, the broadenin g parameter measured in ab- 

which 



Crawford et al. 



1994), 



sorption (b = 2.0 + 0.7 km s" 

on the other hand probes the full column of gas along the 
midplane. 

4. A 2D luminosity profile is computed, where the emission as 
a function of wavelength is derived as a function of cylin- 
drical radius and height above the midplane. The luminosity 
profile Doppler-shifted by the projected velocity is then inte- 
grated through the disk for the line of sight, taking the (also 
Doppler-shifted) line opacity along the line of sight into ac- 
count. 

5. For the line scattering redistribution function we try three 
different cases: 

(a) Complete redistribution, that is, the outgoing photon is 
re-emitted at a wavelength with a probability according 
to the (Doppler-shifted) line profile, independent of the 
wavelength of the absorbed photon. This is a good ap- 
proximation if the Ca n is disturbed (by a collision) be- 
fore emitting again, which is likely not the case in the 
present circumstances because of the strength of the tran- 
sition and the low density of gas. 

(b) Coherent scattering, that is, the re-emitted photon has 
exactly the same wavelength as the absorbed one (only 
Doppler-shifted). This is likely a better approximation, 
but essentially assumes that there is no Doppler-shift of 
the photons due to thermal motion, only the Keplerian 
velocity field. 

(c) Partially coherent redistribution, which assumes co- 
herent scattering in the rest frame of the scatterer. 
Because of thermal motions, this implies a wavelength 
shift between the incoming and outgoing photon on 
top of the Keplerian velocity field. With the thermal 
motions isotropically and Maxwellian-distributed, the 
wavelength shift distribution will depend on the angle be- 
tween the incoming and outgoing direction ( |Unno|1952| l. 
We take this into account by convolving the coherent fine 
luminosity profile with the scattering function, which 
depends on the angle between the star, the volume el- 
ement, and the observer. The thermal motions are as- 
sumed to correspond to the same broadening parameter 
b = 1.5kms~' as before. 

To illustrate the differences between the three assumptions 
we plot a projected height emission profile in Fig. [10] 
Because complete redistribution brings photons in the line 
wings closer to the line centre (on average), it effectively cre- 
ates higher optical depths and strongly suppresses the emerg- 
ing flux from the midplane. For our fiducial model, we used 
partial redistribution. 

The resulting Ca n H profile from the model (with partial re- 
distribution) is shown in Fig. 11 which indeed shows the same 
characteristics as the observations: suppressed emission from the 
midplane and emission visible high above the plane. We empha- 
size that we used precisely the same spatial distribution of gas as 
was inferred from the optically Fei A3860 - we did not attempt 
to carefully fit the Ca ii H emission. Still, the model reproduces 
the observations closely, but not perfectly. Since the model is 
symmetric with respect to the midplane, it does not show the 
strong asymmetries seen in the observations. Furthermore, the 
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Fig. 9: Hydrogen-equivalent gas density inferred from Fei A3860 with parametric fits. Multiply with used abundance fraction to 
obtain the density of any specific element (well-mixed gas is assumed). NE to the left and SW to the right; top panels show cuts 
parallel to the midplane at heights z above the plane and the bottom panels cuts orthogonal to the midplane at cylindrical radii r. 



emission predicted by the model is about a factor 2-3 higher than 
seen in the observations. The discrepancies are not too surpris- 
ing given that we assumed a spatial distribution of gas far beyond 
where it is well constrained by Fe i A3860, e.g. both much higher 
above the midplane and closer in to the star than Fe i is detected. 
Small deviations in the Ca n density can make big differences in 
the Ca ii H line luminosity. This could also explain why the Ca n 
emission extends farther into the SW than the NE, the opposite 
to what is observed for Fe i; it could simply be an optical depth 
effect, through which Ca n is shielded by a higher column den- 
sity towards the star (while the Fe i is optically thin). Overall, 
however, we see no evidence that Fe and Ca are not well mixed; 
the striking difference between the observed intensity maps of 



Fei A3860 and Can H and K is well explained by the difference 
in optical depths between the lines. 

5.2. Comparison between gas and dust profiles 

The general NE/SW asymmetry of the gas disk is, as shown al- 
ready by Brandeker et al. (2004), similar to that of the dust. Our 
parametric model of the gas number density profile can be com- 



pared to the parametric fits of the dust profile by Chen et al 
( 2007]) and |Ahmic et al.| ( |2009| ) to assess how the gas-to-dust 



12 1. Outside 



ratio varies as a function of disk radius (see Fig. 
the break radius the distribution of gas follows that of the dust 
model by Ah mic et al.| (|2009 ), while the results from |Chen et al.| 
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Fig. 10: modelled Can H emission profile, using three different 
photon scattering redistribution functions. The figure shows a 
vertical cut at a distance of 120 AU at the SW side. 
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Fig. 1 1 : modelled Ca n H emission profile, assuming the gas dis- 
tribution derived from Fe i A3860. 



00 AU the gas number density increases towards the star, in 
contrast to the dust. Interpreting the connection between the dis- 
tributions of gas and dust requires more detailed modelling of 
gas production and redistribution, which is outside the scope of 
the current paper (see Xie et al.|2012 for a first attempt at relat- 
ing disk properties to physical production / redistribution mech- 
anisms). 



6. Summary 

Our main results are: 

1 . We have produced the first well-sampled images of the /3 Pic 
gas disk and detected spatially resolved emission from Fe i 
and Ca n. 

2. We confirm the NE/SW asymmetry previously found by 
|Brandeker etaL| ( |20oil >. 

3. We find no evidence for flaring, i.e. the disk scale height in- 
creases linearly with midplane distance to the star. 

4. The disk height profile drops slower than expected for a gas 
disk in hydrostatic equilibrium. We speculate that this might 
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Fig. 12: Number density distribution of gas and dust in the fiPic 
disk midplane, normalised to unity at the break radius, and 
shown for mean of parameters derived for NE and SW. We used 
the one-disk parametric fit from Ahmic et al. ( 2009} . 



2007| l indicate a faster decline of the dust number density. Inside & 



be caused by the dynamical interaction between ions and 
charged dust grains. 
5. The Can H- and K spatial emission profiles indicate that Fe 
and Ca are well mixed throughout the disk; the difference 
to the Fei A3860 profile is caused by differences in optical 
depth. 
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Appendix A: Flux calibration and comparison with 
UVES data 



As mentioned in Sect. 2.2 the photometric accuracy of 
FLAMES/GIRAFFE/ARGUS depends to a great extent on the 
spectral response curve derived from observed standard stars. 
Unfortunately, the stellar models to which the standards are 
compared are very coarse in the blue spectral region observed, 
lacking any detailed spectral features. Moreover, the spectra are 
extracted by summing up over all pixels in the ARGUS array 
when running the EsoRex pipeline, adding sky emission and cos- 
mic ray spikes. Owing to this, a reasonable spectral instrument- 
response curve has to be fitted to some well-chosen points in 
the noisy pipeline-derived instrument response before it can be 
applied to correct our observed target spectra. The fitted re- 
sponse curves from individual observing nights are presented 
in Fig. A.l Such seemingly large differences in efficiency lev- 
els between different observations are unlikely, and must, at 
least partly, be attributed to photometric uncertainties in the 
GIRAFFE spectrograph. Observations of the same standard star 
at three different occasions (green, blue, and cyan lines), e.g., 
suggest a spectral instrument-response that varies by as much 
as 30% at the bluest wavelengths. Taking all these factors into 
account, we estimate a photometric accuracy of at best 20%. 

To examine the relative accuracy of the flux calibration in 
more detail, we plotted the Fe i A3860 and the Ca n H- and K sig- 
nal in overlapping inner and outer observed FLAMES regions, 
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Fig. A.l: Fitted spectral instrument-response curves from ob- 
served standard stars, with colours marking observing dates. 



as well as overlapping data obtained with UVES by Brandeker 
|et al.| ( |2004] >. The Fe i emission seems to be roughly consistent at 
higher latitudes in the inner/outer comparison, but there is a clear 
difference close to the disk midplane, where the strong signal ob- 
served in the inner observations appears to be weaker in the outer 
observations. A similar Fei flux deficit of 30-40% in the outer 
observations is evident in a comparison with UVES data, but the 
situation seems to be reversed in the Can H- and K emission, 
where the FLAMES fluxes at 60 AU are ~1.5 times higher than 
corresponding UVES fluxes, while the data at 120AU match 
fairly well (considering estimated errors for the Ca n extraction 
in that region). A cross-check of the flux-calibrated science spec- 
tra with the available HARPS spectrum of /5Pic allowed us to 
confirm that the slope of the fitted response curves is correct to 
within a few percent, therefore any relative differences between 
the strength of the Fei and Can emission compared to UVES 
data is not caused by errors in the fit. Although the error bars 
assigned to UVES data indicate well-determined fluxes, similar 
difficulties in determining absolute photometric fluxes exist for 
that instrument, therefore any calibration of our data with respect 
to UVES data was not performed. Instead, we assumed that pho- 
tometric errors and a higher seeing during the outer observations 
are likely to blame for the observed flux disagreements, and do 
not affect the qualitative results derived from our modelling (in 
which relative inner/outer discrepancies are also effectively av- 
eraged out). 
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R. Nilsson et al.: VLT imaging of the ji Pictoris gas disk 
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Fig. A. 2: Comparison of Fe i emission in overlapping areas of 
inner and outer NE and SW observations. The pixel columns are 
counted from left to right in the four observed regions. 
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